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We consider the Randall-Sundrum brane-world model with bulk-brane energy transfer where 
the Einstein-Hilbert action is modified by curvature correction terms: a four-dimensional scalar 
curvature from induced gravity on the brane, and a five-dimensional Gauss-Bonnet curvature term. 
It is remarkable that these curvature terms will not change the dynamics of the brane universe at 
low energy. Parameterizing the energy transfer and taking the dark radiation term into account, 
we find that the phantom divide of the equation of state of effective dark energy could be crossed, 
without the need of any new dark energy components. Fitting the two most reliable and robust 
SNIa datasets, the 182 Gold dataset and the Supernova Legacy Survey (SNLS), our model indeed 
has a small tendency of phantom divide crossing for the Gold dataset, but not for the SNLS dataset. 
Furthermore, combining the recent detection of the SDSS baryon acoustic oscillations peak (BAO) 
with lower matter density parameter prior, we find that the SNLS dataset also mildly favors phantom 
divide crossing. 

I. INTRODUCTION 

Numerous cosmological observations have confirmed that the universe is undergoing accelerated expansion. This 
phenomenon was not predicted by conventional cosmology governed by general relativity with the known matter 
constituents. To explain the cosmic acceleration, mysterious dark energy was proposed. There are many dark energy 
models, which can be distinguished in the value and variation of the equation of state (EoS) w during the evolution 
of the universe. The cosmological constant is the simplest candidate of dark energy, whose EoS w = — 1 is located at 
a central position among dark energy models. For quintessence Chaplygin gas 0] and holographic dark energy 
models Q, w always stays bigger than -1. The simplest model with w < — 1 is the phantom models £| but which will 
violate not only the Weak Energy Condition but also Null Energy Condition. In general, dark energy can evolve with 
the change of the EoS w > — 1 to w < - 1 (or vice versa). This transition is called "the phantom divide crossing". 

By current estimates, it seems likely that the phantom divide crossing occurs at the recent epoch [f|, though many 
of these estimations are model dependent. In particular, with an advantageous parametrization (the called CPL 
parametrization) of dynamical dark energy [|| , it has been found @, Q that most of the observational probes indeed 
mildly favor dynamical dark energy crossing the phantom divide at z ~ 0.25. Moreover, in the era of structure 
formation, a highly negative w makes negligible the undesirable dark energy to the total energy density. Hence even if 
the observation still admits w = — 1, it is still useful to construct a more general framework that can permit w < — 1 
0. However, it has been proved that the phantom divide crossing of dark energy described by the minimum coupling 
single scalar field with general Lagrangian is either unstable with respect to the cosmological perturbations or realized 
on the trajectories of the measure zero 10]. So most dark energy models with the phantom divide crossing either 
consist of multiple scalar fields with at least one non-canonical phantom component or must recourse to extending 
gravity theory [12|, HH. The former is usually plagued by catastrophic UV instabilities, however it can be regarded as 
an effective field description following from an underlying theory with positive energies [Til ]. An interesting example 
where the UV pathologies are absent in phantom models is given by [15|. The latter, for example the so called 1/R 
gravity [l6| , is severely constrained by solar system test and by cosmological observation even though it is interesting 
theoretically. A novel model about f(R) gravity which pass through all tests is proposed by [l7|. Other theories with 
the phantom divide crossing are interacting holographic dark energy models [18j |. models with interactions between 
dark matter and dark energy [p|, model with a scale field coupled to the Gauss-Bonnet form [201 ]. model with 
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non-linear term of scalar curvature [2l| , and model considering quantum effects [22J . 

In another area of theoretical cosmology, the brane- world scenario has received a lot of interest. The well known 
Randall-Sundrum (RS) brane- world model 0, |24|, inspired by D-brane ideology in string theory, envisages that our 
four-dimensional universe is a 3 dimensional membrane embedded in the five-dimensional bulk. The standard model 
particles are confined in the brane while the gravity is free to propagate into the bulk. This brane approach provides 
a new way of understanding the hierarchy between the four-dimensional Planck scale and the electro-weak scale. Two 
important generalizations of the RS model have been considered recently. The first is inspired by superstring theory 
which suggests the Gauss-Bonnet (GB) curvature corrections as the first and dominant quantum corrections to the 
Einstein-Hilbert action for a ghost-free theory [25[. The combined action in five dimensions gives the most general 
action with second-order field equation, as shown by Lovelock [26]. The DGP model suggests the second curvature 
correction term to RS model, the four-dimensional scalar curvature term. This induce gravity correction term can be 
interpreted as arising from a quantum effect due to the interaction between the bulk gravitons and the matter on the 
brane [27| . It is natural to study the influence of the GB correction to the DGP model. Indeed, in certain realization 
of string theory, the ghost-free GB term in the bulk action may naturally lead to DGP induced gravity term on the 
brane boundary (28J. It has been suggested that the combination of the two curvature corrections may shed some 
light on the singularity problem in the early universe [29 . |30( • 

However, many brane-world models including the RS and GB types, produce only ultra-violet modifications to 
general relativity, with extra-dimensional gravity dominating at high energy. DGP model may manifest nontrivial 
low energy extra-dimensional gravity, which results in a new branch DGP(+) with late-time acceleration [3l| . Its 
generalization with the brane and bulk cosmological constants [l3| allows w < — 1 and the branch DGP(-) which has 
the right RS limit with w > — 1, but none of these has the phantom divide crossing behavior. 

Another robust way that extra-dimension can affect the low energy evolution of the universe is through the coupling 
between the brane and the bulk. This is analogous to the coupled dark energy scenarios [32| . where the late-time 
accelerating cosmological phase is characterized by a frozen ratio of dark matter/dark energy as a result of the 
interaction of the dark matter with other components, such as scalar fields. It has been proved that the transfer of 
energy between the bulk and the brane may result in the RS model with late-time acceleration [33| and even the 
phantom divide crossing, provided that there is bulk matter 35 1 or an additional dark energy on the brane (34[. 
However, the former can not make cosmological model which is independent of bulk dynamics, and the latter can 
not have the phantom divide crossing by using only geometric effect. The brane-bulk energy transfer has also been 
considered in DGP model for the present universe as a global attractor using fix-points theory [36j and in many other 
different setups [37, 38] . The combined curvature effect and bulk contents effect on the RS model have recently been 
studied in [391 ]. 

The aim of the present work is to study the low-energy cosmological behavior based on an extended RS(II) scenario 
[24| by considering brane-bulk energy transfer and two additional curvature corrections. We obtain a closed system 
of three equations which determines the parameters of the desired Friedmann equation. In the low energy region, it 
is remarkable that the curvature corrections terms do not change the dynamics of the system. By parameterizing the 
energy transfer with the scale factor, we are able to exactly solve the Friedmann equation. When including the effect 
of the dark radiation, which is neglected in |34l. l35j . we show that the phantom divide crossing may be achieved in the 
absence of additional dark energy components on the brane. Furthermore, we fit the model to the two most reliable 
and robust SNIa datasets, the new 182 Gold dataset and the first year Supernova Legacy Survey (SNLS) dataset 
[4lj . respectively complemented by the recent SDSS baryon acoustic oscillations peak (BAO) dataset. 

This letter is arranged as follows: In Sec. II, we establish the most general brane world model with curvature 
correction terms and bulk-brane energy transfer to describe the accelerated expansion. We investigate the equation 
of state of effective dark energy and discuss the possibility of phantom divide crossing. In Sec. Ill, we fit the model 
to the data from dark energy observations. In the last section, we conclude with a brief summary. 



II. BRANE WORLD WITH CURVATURE CORRECTIONS AND BRANE-BULK ENERGY TRANSFER 

Let us consider a braneworld model. For convenience and without loss of generality we choose the extra-dimension 
coordinates y such that the brane is located at y — and the bulk has Z2 symmetry under the transformation y — > —y. 
The most general action which incorporates the induced gravity and Gauss-Bonnet corrections is [301 ] 

S = —-J [ d?x</^g{CEH + oi£gb} + -^s [ d^Xy/^ClG, (1) 

^ K 5 J ^ K A Jbrane 

where g (<?) and K5 (K4) are the bulk (brane) metric and bulk (brane) gravitational constant, respectively. Leh = 
R — 2A is the five-dimensional Einstein-Hilbert Lagrangian with negative cosmological constant A < 0. The Gauss- 
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Bonnet curvature correction term Lqb is written as 

C GB = R 2 - 4RabR ab + RabcdR abcd - 

We can define the Gauss-Bonnet coupling a through string energy scale g s as a = . The induced-gravity 
Lagrangian Ljq = R — 2k|A consists of four-dimensional scale curvature R and brane tension A > 0. We can define 
the induced-gravity crossover length scale by r — — |. For convenience, we will choose the unit that K5 = 1 throughout 
this letter. Note that by setting r = a = 0, we can recover the RS model. The RS model with Gauss-Bonnet correction 
and the DGP model correspond to the case with r = and a = 0, respectively. 

By varying the action in Eq. ([T]) with respect to the bulk metric, we obtain the field equation 

Gab + 2olH A b = T A B\ t otai > ( 2 ) 
where Hab is the second order Lovelock tensor [2(| 

Has — RRab — ^R^Rbc — 2R CD Rabcd + R < a DE ' Rbcde ~ -t9ab^-gb- 
The total energy- momentum tensor Tab\ total 1S decomposed into bulk and brane components 

T AB\ total = T AB\ bu i k + T A B\ brane S(y). 

Here, we use the normalized Dirac delta function, 5(y) = y/ g / gS(y) . The bulk component is 

T AB\ bu i k = ~A 5AB + Tab, 
where Tab denotes any possible energy-momentum in the bulk. The brane component is written as 

TAB\brane = ~^9AB ~ rG AB + Tab , 

where Gab arises from the scalar curvature in Eq. ((T|). The energy momentum tensor Tab represents matter on the 
brane with energy density p and constant equation of state parameter w m . 
The five-dimensional line clement in the bulk is given by 

ds 2 = -n 2 {t,y)dt 2 + a 2 {t,y)j i:j dx l dx j +b 2 (t,y)dy 2 (3) 

where 7^ is a 3-dimensional maximally symmetric metric whose spatial curvature is characterized by k = 0, ±1. In 
this letter, we are interest in spatially flat brane k = 0. We choose the coefficients n(t, 0) = 1 so that t is the proper 
time along the brane. For simplicity, we assume that the fifth dimension is static b = and we set 6=1. 

To determine the Friedmann equation on the brane, we impose the junction condition for a braneworld in Gauss- 
Bonnet gravity. For later use, we define 
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where prime and dot denote the derivative with respect to y and t, respectively. The jump of the (00) and (ij) 
components of the field equation Eq. © across the brane gives 

2[1 + la(H 2 + %]^=rH 2 -u, (4) 

Til a! v! a' ■ 2 • 

2^ + 4^ + 8a— 1$ + l6a-t{H 2 + H) = 3r(H 2 + -H) + v, (5) 
n a n a Q 3 

where 2a', — —2a' and 2n', = — 2n' arc the discontinuities of the first derivatives. H = — is the Hubble constant 

+ ~ + ~ ao 

on the brane and $0 = ®{t, 0). We define 

1 



(/9 + A), v = (w m p-X) 
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One can show that Eq. ([4]) can be written in terms of the square of Hubble parameter H 2 

4[1 + ^(H 2 + \)] 2 (H 2 4> ) = [rH 2 - u} 2 . (6) 

Once we know the behaviors of <l>o and p, Eq. ^ can be solved and its solution will give the desired Friedmann 
equation. 

Let us consider p first. Since the left-hand side of Eq. @ is divergence-free, the total energy momentum tensor is 
conserved in the bulk \7 a Tg I. . , = 0. Its zero component is 

T ° + 3^(T ° - I?) - + 3^)T 5 + Tl' -[p+ 3(1 + w m )Hp]8{y) = 0. 

Integrating around y = and using the Zi symmetry (Tq (i, +0) = — Tq (t, —0)), we can determine the evolution of p 

p + 3(l + Wm )ifp = 2T 05 . (7) 

This implies that the energy conservation law on the brane is broken. Another method to obtain the semi-energy 
conservation law Eq. ([7]) is to solve the junction conditions Eq. (]3|) and Eq. ([5]). Their exact solutions are complicated 
but we can make further simplification by assuming a to be small. This is a reasonable assumption since we are 
interested in the effect of GB correction on the late-time evolution of the universe. The solutions of the Eq. and 
Eq. ([5]) up to the first order in a are 

— = i (rH 2 -u) + \ {rH 2 - u) [r 2 H i + u 2 - 2H 2 (6 + ru)] a, (8) 
ao 2 6 

n' 1 / \ 1 

— t- = - I 2u + v + rH 2 + 2rH) + -{r 3 H 6 + 3riJ 4 (-4 + 2ru + rv + 2r 2 H) (9) 
no 2 V / 6 

+ u[u(8u + 3v) + 6(4 + ru)H] - 3H 2 [4v + u(8 + 5ru + 2rv) + 4r(4 + ru)H]}a. 

Note that Eq. ^ is a cubic equation for the discontinuity ^ which has one real solution, the other two are complex. 
Since we require our cosmological equations to have the right a — > limit [44J, we take only the real solution. 
Substituting Eq. (jHJ) and Eq. ([9]) into the (05) component of the field equation Eq. ((2|) 

3(!l^_^)(l + 4a <I>)=To5, (10) 
n a a 

we recover (up to the first order of a) the same semi-conservation law Eq. ([?])• 

Let us consider the function For brane cosmology with Gauss-Bonnet correction, $ is just the first integral 
satisfying the constraint equation 

$ + 2a $ 2 = - + — , 
6 a 1 

where ^ is the dark radiation term [43} ■ We point out an important observation that the above first integral can be 
recovered from a general differential equation with nontrivial T05 and T55 

^ + 4-^ + 4To5- + 4-(r 55 -A)=0, (11) 
a a a 

where 

ip = 6($ + 2a$ 2 ). 
Notice that we only consider the solution having the right a — * limit 

12a 
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The differential equation Eq. (fTTj) is obtained by substituting ^- from Eq. (fTOf into the (55) component of the field 
equation Eq. |(5J) 



a' n la 
an n 2 a 



(l+4a$)-3$ = T 55 -A. 



Note that for vanishing X05 and T55, the constraint equation Eq. (|lip gives a bulk equation. But for generic values 
of Xo5i the constraint equation depends on the discontinuity of the first derivative on the brane. Moreover, we point 
out that the physical meaning of -0 can be understood as the effective bulk cosmological constant. To see this, we 
rewrite Eq. (fTTj) as 

^i+4-^ 1 +4To5-+4-T 55 -0 (13) 
a a a 

where ipi = ip — A takes the role as the correction to the bulk cosmological constant. We can easily see that, for 
Tq 5 = T55 = 0, the only correction to bulk cosmological constant A is the dark radiation term ~ a~ 4 . For a nontrivial 
bulk, the term 4Tos ^- + 4£f T55 in Eq. (fTTj) will give further corrections to A. We will come back to discuss this point 
later. 

Generally, Eq. © has three solutions of H 2 , but only two of them are left for small a: 



H 2 = 



6 + 3ru =F 1 



3r 2 ^ 9r 4^ 



{32[=F9r 2 it + 6ru(/> qp 2<Xt6 + 4>)} ± 8r 2 (±18 + 3™ + <j>)ij) ± rV 2 }, (14) 



where we have used Eq. (fT2")l and = vGy 6 + 6ru — r 2 ip. These two branches recover the two branches of DGP 
model when a — > 0. One can find three equations (JT]), (JTTJ) and (H3J) consisting of a closed system for variables p, 
$0 and H, which determines the evolution of the universe (provided the bulk energy- momentum tensor is known). 
However, these equations are entangled in a complicated way. In the following, we will seek the simplest low-energy 
effective theory which can be analytically solved. 

We would like to make some assumptions here. First, we consider constraints on p and H 2 . In low energy region, 
we can assume 

A A , A . . 

p < A, /?<-,/?<—, I/ 2 <-. (15 

This low energy region is self-consistent and can be easily realized. The simplest example is the case where the brane 
tension and the bulk cosmological constant are very large. Moreover, it should be noticed that in the a — > and 
r — > limits Eq. (|15p is consistent with RS low energy region 

p< A, H < A, \p~H 2 , A- A 2 (16) 

where the last constraint comes from the well known RS fine-tuning mechanics. Noticing that the present Hubble rate 
Hq is much smaller than the string energy scale g s and H 2 /Hq is close to 1 at late time, we have aH 2 <C 1. Hence 
the region p <C -tj does not suggest new energy constraint beyond Eq. (|16p . Similar consideration can be applied to 
the rH 2 <C A constraint in Eq. (fT5"j) if we set rH ~ 1 in DGP(+) brane [29(. Second, to derive a cosmological system 
that is largely independent of the bulk dynamics, usually one can assume that on the brane the contribution of T55 
relative to the bulk vacuum energy is much less important than the brane matter relative to the brane vacuum energy, 
or schematically j. Thus, T55 can be omitted in the low energy region. Third, the semi-conservation law Eq. 

(JT]) implies T05 ^ Hp in general. Fourth, at late time, the dark radiation may be assumed to be much smaller than 
bulk cosmological constant. By taking in to account the above assumptions, we find that Eq. (|13p can be simplified 
as 

^1 + 4ffVi + 4T 05J 4 - 0, (17) 

where A is a constant defined by 



It should be noticed that Tq$A <C HA. Returning to the discussion about the correction to the bulk cosmological 
constant. We can now conclude that the correction to the bulk cosmological constant ipi is much smaller than the 
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bulk cosmological constant, ■01 A. This is reminiscent of p <C A on the brane, and can be understood as bulk-brane 
duality. Thus, Eq. (|14[) can be expanded in terms of p and ipi (keeping only the first order terms): 

H 2 = Bp + Cipi + D, (18) 

where we absorb A and A into the constant coefficients 

B = iSr + ^ E T ~ 9r [ ±8rA ( 8 + 3rA ) ± 32A ( 4 + rA ) ± 3^A 2 ]}, 

C = ±- + -^[-96(£A qp 6) - 32rA(£ t9) + 8r 2 A(2^ =p 27) T 3r 4 A 2 ], 
hi E a r z 

D = 6 ^ E + rX + -^-{64Er\ T 64£(6 T E) T 32r 2 A 2 ± 8r 2 [ T A(-18 T E + rA) ± r 4 A 2 ]}, (19) 
3r z 9£r 4 



with J5 = ^6(6 + 2rA - r 2 A). 

Observing the three equations Eq. 0, Eq. (|T7|) . and Eq. (|T8|) . we find that the curvature corrections only 
affect the constants A, B, C, D, E, and they will not vanish in the a — > and r — > limits. This remarkable result 
immediately tells us that the curvature corrections will not change the form of the desired Friedmann equation, and 
the corresponding evolution of the universe. 

To obtain the explicit solution we need to know the form of energy transfer Xo5 ■ Unfortunately, it is not yet available 
and obviously depends on the mechanism which produces the energy transfer. We consider the ansatz T05 = THa" 
(T is a constant) which was used in [34], HH for RS brane world. A justifcation for the ansatz has been analyzed in 
[33 ] through a simple model where the bulk content is a relativistic fluid, slowly moving along the fifth dimension. 
Recently, the ansatz has also been used in j4^|. Then we find 

IT 

p = Fa' 3 + a" (20) 

3 + v 

where we take w m = for dark matter, and 

where G is an integration constant. Notice that the last term in Eq. (|20|) denotes the energy flow into (j^j > 0) or 
out of (^2_ < 0) the brane. Substituting Eq. $IU§ and Eq. into Eq. US]), we have 

H 2 = Qom,a~ 3 + ftova v + fWa~ 4 + D, 

where 

^om = BF, Clov — B— h C — —j- — - — -, Clod = CG. 

6 + v 81(4 + v) 

To impose the vanishing effective cosmological constant on the brane, we choose D = 0. It just recovers RS fine-tuning 
when a — ► 0. Finally, we can write 

H 2 = n 0m a" 3 + n Qd a- 4 + fl 0v a v . (22) 

This is similar to the well known result in RS model where the Friedmann equation can be generalized by adding the 
term depending on the brane-bulk energy flow. The curvature corrections only affect the coefficients n 0( „, fio<i and 
flov However, it should be pointed out that, they may play an important role, for example, when the constant G 
in dark radiation term and the parameter T which characterizes the brane-bulk energy exchange are not big enough 
to provide desired nontrivial cosmological behavior. In particular, the dark radiation term is usually neglected if one 
considers the large scale factor at late time. In RS model, the reason that we can keep it is that indeed we do not 
know the constant G which reflects the bulk geometry (Notice that it takes the role of the bulk black hole mass in 
a Schwarzschild-AdSs geometry 44]). In our model, furthermore, the curvature corrections embodied in constant C 
may strengthen the need for keeping the dark radiation term. For an explicit example, we assume G <SC 6FXa so that 
the dark radiation term can be omitted f2odfl~ 4 *C fiom& 3 hi RS case where C = g, B = A. For the dark radiation 
term to be important fiodO -4 ~ fiomO -3 , we need G ~ <C 6F\a i.e. |j <C 6A. Considering the induced gravity 
correction and using lower branch for simplicity, we find that the condition is satisfied when r 3> 54 3 ^ A , which does 
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FIG. 1: The EoS w e ff (red line) with prior zt = 0.2, v = 1 and deceleration factor q (green line) with prior zt = 0.2, v = 1, 
Hom = 0.25 versus redshift 2. One can find that w e ff crosses -1, and the q crosses at z ~ 0.5 and the order of magnitude 



not violate the low energy region Eq. (fT5|) . Another example can be given by considering the relationship between the 
accelerated expansion and the brane-bulk energy flow. The accelerated expansion is characterized by the deceleration 
parameter q — —jjt- Using a — and taking iJ 2 | z _ Q = 1, we can rewrite the Friedmann equation as 

h 2 = Q m(i + zf + n 0v (i + zy v + (1 - n Qm - o 0l ,)(i + zf 

where we have absorbed into f^om and flo v . The deceleration parameter can be given as 

1 dlogH 2 
q=-l+-- 



2 dlog(l + 2) 



Now we will omit the terms which decrease more quickly than matter density at late time, then the dark radiation 
term is absent and one can find that the present accelerated expansion q\ z=0 < needs v > — 3 and tto v > > 0. 
This suggests that in RS case (noticing A > in RS(II) brane world), the bulk energy must flow into the brane (the 
later term in Eq. (|20[) is positive, i.e. T > 0). However, if we consider the curvature correction (still the induced 
gravity correction and using the lower branch), we find that accelerated expansion may be achieved when the bulk 
energy flows out of the brane, which only needs |||| < gz^n ^ that can be easily realized. 

Let us consider whether or not this model is permitted to cross the phantom divide and what is favored by fitting 
the model to observational datasets. We find that the combination of the last two terms in Eq. (|2"2"]) may achieve the 
phantom divide crossing, without the need of any other dark energy components. The EoS of effective dark energy 
can be written as [45( 

1 , 1 d\og[H 2 - n 0m (l + z) 3 } 

W eff = -1 



3 d\og(l + z) 



It is easy to see that w e f / is a constant if the dark radiation term is omitted Qq v = 1 — Slo m or the brane-bulk energy 
transfer is trivial f2ot> — 0. But if both of them are taken into account, we may achieve the phantom divide crossing 
in very broad parameter space. Explicitly, let us assume w e ff = — 1 at zt, then 

™ e //U =0 - 1+ 3b + 4(l + z T )4+-]- 

When zt > 0, one can find w e ff\ z=0 < —1 if v > 0. This implies that the brane-bulk energy flow increases with the 
expansion of the universe. For an explicit example of the phantom divide crossing and current accelerated expansion, 
see Fig. 



III. OBSERVATION OF THE PHANTOM DIVIDE CROSSING 



Now we use the observational datasets to test our cosmological model. We will use the new 182 gold supernova la 
data (0 < z < 1.76) and the first year Supernova Legacy Survey (SNLS) dataset (0 < z < 1), combined with the recent 
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BAO measurement from SDSS to fit our model. Recently these data have been widely used in the fittings of different 
braneworld cosmological models |47|. There are other dark energy observational probes, including the 3-year WMAP 
CMB shift parameter [48[ , the X-ray gas mass fraction in clusters [HI and the linear growth rate of perturbations at 
z = 0.15 as obtained from the 2dF galaxy redshift survey [50]. However, since our model is effective at low energy, 
we will not use the 3-year WMAP CMB shift parameter focusing on the high redshift region. Besides, for simplicity, 
we do not adopt other probes of dark energy which have large relative errors compared with Snla, SNLS, CMB and 
BAO probes [|. 

For the supernova la data, the measured quantity is the bolometric magnitude m 

m = M + 51og 10 (Z3i) 
where M is the Hubble-parameter-free absolute magnitude 



M is the absolute magnitude, and 



M = M-51og 10 (^) + 25. 



D L = (l + z) I dz' 



(i 



H(n 0m ,$i 0v ,v) 



is the Hubble free luminosity distance (HqcIl/c). The data points of the Gold dataset are given after implementing 
correction for galactic extinction, if -correction and light curve width-luminosity correction, in terms of the distance 
modulus 

uC° ld (z t ) = m% s ld ( Zi ) - M. 

For SNLS datasets, also presents for each point, the stretch factor s used to calibrate the absolute magnitude and the 
rest frame color parameter c which mainly measures host galaxy extinction by dust. Thus, the distance modulus in 
this case depends apart from the absolute magnitude M, on two additional parameters a and (3 

u s Z LS { Zl ) = m s o a LS { Zl ) - M + a{ Sl - 1) - p Cl . 

Let us define the theoretical distance modulus 

uth(zi) = m th {zi) - M = h\og w {D L ) - 51og 10 (— p— ) + 25. 

We shall assume that the supernova la measurements come with uncorrelated Gaussian errors of (including flux 
uncertainties, intrinsic dispersion of Snla absolute magnitude and peculiar velocity dispersion) in which case the 
likelihood function is given by the % 2 distribution 

2/ o ^ [Uobs(Zj) ~ U th (Zj)} 2 
X ("Orrn "0t» = > , 2 ( 23 ) 

where N = 182 for Gold datasets, and N = 115 for SNLS. 

For BAO measurement, we shall use the model independent measurement of the parameter [5l| 



O 2 

"Om 



Ho 
H 



-dz 



0.35 J H 

to construct an additional term in the % 2 equation Eq. (|23p 



0.469 ±0.017 



2 _ [A(n Om ,n Ov ,v)-0A69} 2 
Xbao 0.017 2 

The theoretical model parameters are determined by minimizing the % 2 and \ 2 + Xbao- 

In order to investigate the dependence of the resulting best fits on the prior of Slo™, we will consider two cases fiom 
= 0.2 and = 0.3 instead of marginalizing over f2om- The range between the two cases includes the current best fit 
value of fi 0m based on WMAP and SDSS which is f2 0m = 0.24±0.02 J52|. Considering the errors using the covariance 
matrix method (53|, we show the best fit form of w e ff for each dataset category in Fig. @. The corresponding \ 2 
contours in the remaining two parameter space is shown in Fig. ([3]). 
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FIG. 2: The best fit form of w e ff(z) for different dataset category for both Slom = 0.2 and fiom = 0.3. The categories are: 
Goldl82 dataset and it with BAO (row 1), SNLS and it with BAO (row 2). The dashed line in each panels represents the 
phantom divide. 
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FIG. 3: The 68% and 95% confidence contours in the Slov — v parameter space for each dataset category for both Q,Q m = 0.2 
and !!o m = 0.3. The cross point of dashed lines in each panels represents ACDM. 

IV. SUMMARY 

We have studied the cosmological dynamics of the general RS brane world scenario with brane-bulk energy transfer 
and two curvature correction terms. For small Gauss-Bonnet coupling, we obtain a closed system of three equations 
having two branches that correspond to the two branches in the DGP model. They describe the evolution of the 
Hubble rate, the energy density and the time-dependent effective bulk cosmological "constant". We find that in the 
low energy region Eq. (fT5|) , these two branches have equivalent dynamics as that of the RS model with the brane-bulk 
energy exchange. Furthermore, we have shown that the phantom divide crossing can be achieved when the dark 
radiation term is presented and the brane-bulk energy flow increases with the expansion of the universe, without 
the need of any additional dark energy components on the brane [34j . or bulk matter [3o| . or even exotic phantom 
material From the fitting, it has been revealed that the model indeed has a small tendency of phantom divide 
crossing for the Gold dataset, but not for the SNLS dataset. This is consistent with the analysis in 0, [f| for CPL 
parametrization of dynamical dark energy. We have further shown that the BAO constraint with the lower matter 
density prior mildly changes the tendency of SNLS dataset and favors the w crossing —1. 
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